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Dredging  Research 

T echnical  Notes 


Tidal  Constituent  Database  —  East  Coast, 
Gulf  of  Mexico,  and  Caribbean  Sea 
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This  technical  note  describes  a  database  of  tidal  elevation  boundary  con¬ 
ditions  generated  in  support  of  "Long-Term  Fate  of  Dredged  Material  Dis¬ 
posed  in  Open  Water"  research  of  the  Dredging  Research  Program  (DRP), 
being  conducted  by  the  U.S.  Army  Engineer  Waterways  Experiment  Sta¬ 
tion  (WES).  The  database  described  allows  the  user  to  manually  generate 
time  series  of  tidal  elevations  or  to  use  a  program  to  access  the  full 
database  to  generate  time  series  of  both  tide  elevations  and  currents  for 
any  location  along  the  U.S.  east  coast.  Gulf  of  Mexico,  and  Caribbean  Sea. 
Although  the  capability  to  generate  these  time  series  was  developed  to  pro¬ 
vide  input  to  the  long-term  fate  and  stability  model  LTFATE,  the  gener¬ 
ated  time  series  can  be  used  for  any  application  requiring  tidal  forcing 
data. 

Background 

The  long-term  fate  research  has  been  concerned  with  developing  tech¬ 
niques  to  predict  the  long-term  fate  of  dredged  material  after  it  has  been 
deposited  in  open  water  on  the  ocean  floor,  that  is,  to  address  the  question 
whether  a  dredged  material  disposal  site,  either  existing  or  proposed,  is 
dispersive  or  nondispersive  (Scheffner  1992).  If  the  site  is  dispersive,  an 
additional  capability  of  the  model  is  to  estimate  the  rate  of  erosion  and 
fate  of  the  material.  Because  sediment  is  primarily  eroded  and  transported 
as  a  function  of  waves  and  currents,  the  approach  is  to  construct  databases 
of  site-specific  information  that  can  be  used  as  input  to  coupled  hydrody¬ 
namic,  sediment  transport,  and  bathymetry  change  models  for  predicting 
the  long-term  behavior  of  disposal  sites.  In  the  DRP,  attention  has  been  fo- 
'  cused  on  the  development  of  the  wave,  tidal,  and  storm  surge  components. 


The  wave  component  of  the  database  provides  the  capability  for  generat¬ 
ing  time  series  of  wave  height,  period,  and  direction  for  any  location  at 
which  a  WES  Wave  Information  Study  hindcast  is  available.  The  wave 
simulation  capability  is  described  in  Borgman  and  Scheffner  (1991).  A 
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database  of  regional  tropical  and  extratropical  storm  surge  elevation  (and 
current  hydrographs  and  stage-frequency  relationships)  is  currently  under 
development.  This  technical  note  describes  the  tidal  component. 

Additional  Information 


For  additional  information,  contact  the  author.  Dr.  Norman  W.  Scherr- 
r.er,  (601)  634-3220,  or  Mr.  E.  Clark  McNair,  Jr.,  manager  of  the  Dredging 
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The  tidal  contribution  to  the  overall  database  has  been  completed  for 
the  entire  east  coast  of  the  United  States,  Gulf  of  Mexico,  and  Caribbean 
Sea  (Westerink,  Luettich,  and  Scheffner  1993).  The  tidal  component  is 
both  co-tidal  chart-based  for  manual  simulation  of  time  series  and  PC- 
based.  Both  versions  are  described  briefly  in  this  technical  note. 


This  tidal  database  consists  of  precomputed  amplitudes  and  Greenwich 
epochs  corresponding  to  the  eight  primary  aiumai  and  semidiurnal  astro¬ 
nomical  constituents  (Kp  Oj,  Pj,  Qj,  M2,  S2,  N2,  and  K2).  The  constituents 
are  based  on  a  6-month  simulated  tidal  time  series  computed  by  the  long¬ 
wave  hydrodynamic  finite  element  model  ADCIRC  (Luettich,  Westerink, 
and  Scheffner  1992).  This  model  was  developed  for  the  DRP  for  the  spe¬ 
cific  application  to  long-term  computations  over  very  large  computational 
domains. 


The  tidal  constituents  were  computed  by  an  eight-constituent  harmonic 
analysis  of  model  elevation  and  current  time  series  at  each  node  of  the 
computational  grid  shown  in  Figure  1.  Boundary  conditions  for  the  simu¬ 
lation  were  provided  at  the  mid-Atlantic  Ocean  boundary  according  to 
Schwiderski's  (1980)  global  ocean  database.  Use  of  the  DRP  tidal  cfatabase 
will  allow  the  user  to  generate  time  series  of  tidal  elevations  and  currents 
at  any  location  within  the  computajtional  domain  for  any  time  period — 
past,  present,  or  future.  Although  the  intent  of  the  simulation  capability 
is  to  provide  time  series  input  to  the  long-term  fate  model  LTFATE,  the 
generated  data  can  be  used  for  any  application  requiring  tidal  forcing 
(tidal  circulation  studies). 


The  following  section  desciibes  the  procedure  used  to  generate  time  se- 
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and  epochs.  Then,  descriptions  are  given  of  applying  the  procedure  with 
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elevation  time  series  or  to  use  a  computer-based  program  to  generate  sim¬ 
ulated  forecast  or  hindcast  time  series  of  surface  elevations  and  currents. 
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Figure  1.  Tidal  database  computational  domain 


Harmonic  Reconstruction  of  Tides 

The  tidal  elevation  and  current  time  series  at  any  location  can  be  writ¬ 
ten  as  a  function  of  known  harmonic  constituents  according  to  the  follow- 
-  ing  general  relationship: 
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total  number  of  tidal  constituents  in  the  series  reconstruction 

node  factor  for  reducing  mean  amplitude 

mean  amplitude  of  tidal  constituent  n 

speed  of  constituent  n  in  degrees/unit  time 

time  measured  from  some  initial  epoch 

value  of  the  equilibrium  argument  for  constituent  n  at  some 
location  when  t  =  0 

local  epoch  of  constituent  n 


In  the  above  formula,  the  tide  is  represented  as  the  sum  of  a  coefficient 
multiplied  by  the  cosine  of  its  respective  arguments.  A  finite  number  of 
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for  the  arguments  are  generally  computed  through  least  squares  analyses 
of  prototype  data  (Dronkers  1964).  The  National  Oceanic  and  Atmo¬ 
spheric  Administration's  National  Ocean  Survey  .(NOS)  generally  uses  37 
constituents  in  its  published  harmonic  analyses.  These  results  are  usually 
based  on  an  analysis  of  a  minimum  of  1  year  of  prototype  data.  As  an  ex¬ 
ample,  a  listing  of  the  NOS  constituents  and  the  corresponding  ampli¬ 
tudes  and  local  epochs  for  the  tide  station  located  at  Duck,  NC,  is  shown 
in  Table  1.  A  thorough  description  of  each  constituent  is  provided  in 
Dronkers  (1964)  and  Schureman  (1958). 


Most  constituents  in  Table  1  are  associated  with  a  subscript  that  indi¬ 
cates  the  approximate  number  of  cycles  per  solar  day  (24  hr).  The  constit¬ 
uents  with  subscripts  of  2  are  classified  as  semidiurnal  constituents  and 
produce  a  tidal  contribution  that  occurs  approximately  twice  a  day.  Diur¬ 
nal  constituents  occur  approximately  once  a  day  and  have  a  subscript  of 
1.  Symbols  with  no  subscript  are  termed  long-period  constituents  and 
have  periods  greater  than  a  day. 

The  majority  of  constituents  shown  in  Table  1  can  be  neglected  along 
the  east  coast.  Gulf  of  Mexico,  and  Caribbean  Sea.  For  example,  in  the 
modeled  domain  of  Figure  1,  over  90  percent  of  the  tidal  energy  can  be 
represented  by  the  M2,  S2,  N2,  and  K2  semidiurnal  and  Kj,  Oj,  P-j,  and  Qj 
diurnal  constituents.  For  example,  these  eight  constituents  account  for 
94  percent  of  the  energy  of  a  signal  computed  for  Duck,  NC,  with  the  full 
37  NOS  constituents.  For  this  reason,  the.  above  eight  constituents  are 
used  to  define  the  tidal  contribution  to  the  east  coast.  Gulf  of  Mexico,  and 
Caribbean  Sea  database.  In  other  locations  in  the  world,  many  more  tidal 
constituents  may  be  needed  to  adequately  represent  the  tide.  For 
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Table  1.  NOS  Tidal  Constituents  for  Duck,  North  Carolina  j 
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205.48  | 
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0.0140 
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0.0120 
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0.0170 
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example,  to  define  the  tide  at  Anchorage,  AK,  over  100  constituents  are 
necessary. 


Two  categories  of  tidal  constituents  are  computed — those  that  represent 
the  elevation  of  the  water  surface  and  those  that  specify  time.  For  exam¬ 
ple,  the  value  for  Hn  in  Equation  1  is  the  constituent  amplitude  and  is  a 
function  of  both  location  and.  variations  arising  from  changes  in  the  lati¬ 
tude  of  the  moon's  node.  The  moon's  nodal  effect  is  reflected  by  the  intro¬ 
duction  of  the  node  factor  /„,  which  modifies  each  constituent  amplitude 
to  correspond  to  a  specific  time  period.  Midyear  values  are  usually  speci¬ 
fied  for  a  given  time  series  reconstruction  because  node  factors  vary 
slowly  in  time.  Midyear  values  for  each  constituent  of  Table  1  are  pub¬ 
lished  for  the  years  of  1850  to  1999  (Schureman  1958). 


The  second  category  of  arguments  specify  the  timing  of  the  individual 
constituent  high-water  mark  with  respect  to  both  local  time  and  global 
time.  These  arguments  are  based  on  the  fact  that  phases  of  the  constitu¬ 
ents  of  the  observed  tide  do  not  coincide  with  the  phases  of  the 
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corresponding  constituents  of  the  equilibrium  tide.  For  example,  a  high 
tide  does  not  occur  directiy  beneath  the  moon.  There  is  a  lag  between  the 
high-water  phase  of  the  argument  (that  is,  location  of  the  moon)  and  the 
observed  time  of  high  water.  This  lag  is  referred  to  as  the  epoch  of  the 
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ing  high  tide  is  shown  in  the  schematic  diagram  of  Figure  2.  In  this  fig¬ 
ure,  time  is  increasing  to  the  right,  and  k  represents  the  phase  lag  or  time 
required  for  the  water  surface  to  reach  high  water  (HW)  following  the 
crossing  of  the  moon  (M).  Because  the  water  does  not  respond  exactly  ac¬ 
cording  to  theory,  the  value  of  k  is  computed  as  the  sum  of  the  theoretical 
argument  (Vq  +  u )  and  the  actual  observed  phase  angle  £  at  some  time  T. 

Because  constituents  can  be  considered  harmonic  (can  be  expressed  as  a 
cosine  function  whose  argument  increases  linearly  with  time),  the  value  of 
k  is  relatively  constant  at  every  location.  That  is,  the  value  of  K  represents 
the  actual  iag  between  the  tidal  potential  and  the  following  high  tide  as  a 
function  of  observational  data  adjusted  to  reflect  time  to  equilibrium  the¬ 
ory.  Therefore,  the  value  of  k  can  be  computed  from  the  value  of  £  de¬ 
rived  from  prototype  data  measured  at  any  time  and  the  corresponding 
adjustment  according  to  (Vq  +  u)  for  that  location  at  that  time.  Values  of 
the  equilibrium  argument  for  the  constituents  of  Table  1  relative  to  the 
passing  of  the  tidal  potential  at  the  Greenwich  meridian  are  published  for 
each  calendar  year  from  1850  through  2000  (Schu reman  1958). 


Phases  of  a  given  tidal  constituent  in  different  parts  of  the  world  are 
not  directly  comparable  with  respect  to  the  local  epoch  k  because  K  is  a 
function  of  the  longitude  of  the  specific  location.  However,  an  adjusted 


Figure  2.  Tidal  phase  relationships 
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epoch  can  be  computed  which  is  independent  of  both  longitude  and  the 
time  meridian.  This  epoch  has  been  designated  as  the  Greenwich  epoch, 
G,  and  is  related  to  the  local  epoch  k  of  Equation  1  according  to  the 
following: 


Greenwich  q)och  (G)  ~  k  +  pL  (2) 

where  p  is  the  coefficient  indicating  the  number  of  cycles  per  day  (1  for  diur¬ 
nal  and  2  for  semidiurnal)  and  L  is  the  longitude  of  the  station.  The  epochs  of 
the  DRP  database  are  Greenwich  epochs  because  the  time  series  generated  by 
the  ADCIRC  model  are  computed  with  boundary  conditions  referenced  to 
the  Greenwich  epoch  (Schwiderski  1980). 

The  following  section  demonstrates  use  of  the  DRP  database  to  man¬ 
ually  generate  an  equilibrium  surface  elevation  time  series  of  data  at  any 
location  of  the  computational  domain  shown  in  Figure  1.  This  application 
provides  the  user  a  rapid  capability  of  estimating  tidal  elevation  data  at 
any  location  without  use  of  computer  support. 


Tidal  Surface  Elevation  Time  Series — Manual 
Reconstruction 

One  application  of  the  database  is  to  provide  a  means  of  specifying  a  re¬ 
alistic  tidal  elevation  at  some  specific  location  but  not  to  provide  a  tide 
prediction  that  is  accurate  in  both  magnitude  and  time  of  high  and  low 
tide.  Because  the  precise  time  of  arrival  of  the  ride  is  not  important  in  a 
long-term  disposal  site  stability  application,  only  the  Greenwich  epochs 
are  used.  Therefore,  the  following  example  demonstrates  the  generation 
of  an  equilibrium  tide  in  which  the  nodal  factor  /  is  1.0,  the  equilibrium  ar¬ 
gument  (Vq  +  u )  is  0.0,  and  the  local  epoch  k  is  replaced  by  the  Green¬ 
wich  epoch  G.  However,  if  a  tidal  hindcast  or  prediction  is  desired  and 
computer  resources  are  not  available,  guidelines  for  computing  the  node 
factor  and  equilibrium  arguments  are  given  in  Schureman  (1958),  and  the 
values  of  k  can  be  computed  from  the  values  of  G  according  to  Equa¬ 
tion  2.  If  computer  resources  are  available,  guidelines  for  generating  time 
series  of  tidal  elevations  and  currents  are  provided  in  the  next  section. 

To  reconstruct  the  tide  for  any  location,  the  values  of  the  amplitudes 
and  Greenwich  epochs  for  a  particular  location  must  be  extracted  from  the 
database.  For  the  manual  reconstruction  approach,  detailed  co-tidal  charts 
are  provided  in  Westerink,  Luettich,  and  Scheffner  (1993).  An  example  of 
the  M2  charts  for  the  cast  coast  is  shown  in  Figure  3  for  amplitude  and  in 
Figure  4  for  Greenwich  epoch.  The  steps  described  in  the  following  para¬ 
graphs  are  performed  to  generate,  or  resynthesize,  a  tidal  signal. 
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Step  1 

Interpolate  amplitudes  and  phases  for  the  eight  astronomical  constitu¬ 
ents  from  the  co-tidal  charts.  For  the  Duck,  NC,  example  presented 
above,  the  amplitudes  and  Greenwich  epochs  shown  in  Table  2  were  ex¬ 
tracted  from  Westerink,  Luettich,  and  Scheffner  (1993).  The  constituent 
speeds  shown  in  the  table  are  readily  available  from  sources  such  as 
Schureman  (1958)  or  Dronkers  (1964). 


|  Table  2.  Harmonic  Arguments  for  Duck,  North  Carolina 

Constituent 

Amplitude  H,  m 

Epoch  G,  deg 

Speed  to,  deg/hr 

Kt 

0.092 

179 

15.0410686 

Oj 

183 

13.9430356 

Pi 

175  14.9589314 

Qi 

179  :  13.3986609  j 

n2 

0.14 

354 

28.4397296  ! 

m2 

0.63 

14 

28.9841042  [ 

S2 

0.11 

5  7 

30.0000000  ■ 

k2 

0.037 

48 

30.0821372  j 

Step  2 

Compute  the  tide  signal  according  to  Equation  1.  The  resynthesized 
tidal  elevation  signal  for  the  data  between  days  50  and  60  is  shown  in  Fig¬ 
ure  5. 

Tidal  Surface  Elevation  and  Current  Time 
Series — Computer  Generation 


Tidal  time  series  based  on  the  surface  elevation  and  current  (u  and  v 
components)  amplitude  and  the  Greenwich  epoch  database  arc  available 
for  mainframe  or  PC  applications.  The  computer  program  that  accesses 
these  data  has  the  capability  of  computing  equilibrium  arguments  (Equa¬ 
tion  1)  so  that  hindcast  or  forecast  time  series  can  be  generated.  There¬ 
fore,  either  equilibrium  time  series  or  time-referenced  time  series  can  be 
generated.  Input  to  the  program  includes  equilibrium  or  hindcast /fore¬ 
cast  option,  length  of  desired  time  series  and  time  increment  between  data 
-  points,  location  latitude  and  longitude  if  a  hindcast /prediction  is  desired, 
and  starting  time,  in  hour,  day,  month,  and  year. 

The  program-generated  prompts  for  a  30-day  tidal  hindcast  at  Duck, 

NC  (latitude  36.1833,  longitude  75.7333),  beginning  at  0000  hr  on  1  July 
1993,  are  shown  in  Figure  6.  The  generated  time  series  is  shown  in  Fig¬ 
ure  7.  Node  factors,  equilibrium  argument  amplitudes,  and  epochs  (both 
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Figure  7,  Computer-generated  tidal  elevation  and  current  time  series  for  Duck,  NC, 

July  1993 
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k  and  G),  and  time  series  used  to  generate  Figure  7  are  contained  in  a  file 
named  tide.plot,  shown  in  part  in  Figure  8. 
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A  database  of  tidal  elevation  and  current  astronomical  arguments  has 
been  completed  for  the  east  coast  of  the  United  Staies,  Guif  of  Mexico, 
and  Caribbean  Sea.  Co-tidal  charts  of  surface  elevation  amplitudes  ana 
Greenwich  epochs  for  the  east  coast  and  Gulf  of  Mexico  are  provided  m 
Wes  ter  ink,  Luettich,  and  Scheffner  (1993).  That  publication  provides  the 
capability  of  generating  tidal  elevation  time  series  for  any  location  in- 
eluded  in  the  charts  &nd  docs  not  require  the  use  of  a  computer.  The  full 
database  of  surface  elevation  and  current  amplitudes  and  epochs,  as  well 


h2crpRsO:larry$  pq  tide.plot 
1,0  HOUR  TINE  INCREMENT  1 


75,73330  36.18330 


30  MY  TIDE  PREDICTION  STARTIKfi:  HR- 0.00,  DAY-  1,  MONTH-  7  YEAR- 1333 

CONST  NODE  EQ  ARB  ELECTION  EAST  VELOCITY  NORTH  VELOCITY 
NAME  FACTOR  (DEG)  <MP(M)  G(DEG)K(DEG)  AHP(H)  6(DE6)K(DES)  AMP(H)  G(DEfi)K(DES) 


X 1 
01 
FI 


H2 

S2 

K2 


0.37400 
0.95771 
i  .00000 
8,95771 
1.01231 
1.01231 
1.00000 
0 .51871 


137,79 

242.59 

178.34 

183.51 

25.46 

Ml  V*J 

w.ys 

360.00 

215.28 


0,0923 

0.0659 

0.0381 

Hr01?l 

0.133? 

8.6293 

0,1119 

0.0366 


173.4 
183.6 

175.5 
173.3 
354.0' 

13.7 

57.2 

48.4 


103.6 

107.8 

35.8 

103.5 

202.6 

ttt.  2 

265,7 

256.9 


0,0025 

0.0017 

0.0010 

n.fifio? 

0.0071 

u.utfw 

0.0063 

0.0020 


327.5 

338.1 

312.2 

340.6 

114.7 

132.8 

173.3 

163.4 


251.8 

262.4 

236.5 

264.8 

323.3 

341.3 

21.3 
11.9 


0,0048 

0.0033 

0.0020 

n.nniM 

0.di73 

3.0743 

0.0161 

0.0052 


141,2 

154.0 

125.4 

154.8 

282.7 

299.7 
342.6 

332.8 


65.5 
78.2 

49.6 
73.(1 

131.3 

148.3 
191.1 

181.3 


TIKE  (HRS)  AMPLITUDE  (N)  EAST  VELOCITY  (M/SEC)  NORTH  VELOCITY  (H/SEC) 


0.000 

-0,66711 

0.00847 

-0.00119 

1.000 

-0.50549 

-0.00623 

0.03571 

2.000 

-0.23206 

-0.01866 

0,06231 

3,008 

0.07807 

-0.02563 

0,07179 

4.0DB 

0.34191 

-0,02533 

0.06172 

5.000 

0.4B34? 

-0.01782 

0.03456 

6.000 

0.48150 

-0.00500 

-0.00282 

7.060 

0.31346 

0.00584 

-0.04O92 

8.000 

0,04525 

0,02287 

9.000 

-0.26915 

0.03063 

-0.0’8220 

10.000 

-0.53979 

0.03054 

-0.07415 

11,986 

-0,63186 

n_A?24* 

-0,94733 

12.000 

-0.67884 

0.00978 

-0.007a) 

13.000 

-0.00694 

14,908 

-0,17586 

-0.02272 

0,07078 

15.000 

0.20672 

-0,03378 

0.03081 

i  r  aam 

lo.vou 

0.565GG 

r>  in 

-U.UJ/Hi 

A  AJVSt  -1 

0.U7U1/ 

17.000 

0.81500 

-0,03272 

0,06891 

18.000 

0.89741 

-0.02077 

0.03211 

19.008 

0.79224 

-0.00440 

-2.811SS 

20.000 

0.52379 

0.0124B 

-0.05111 

2f.doo 

0.15483 

0.02583 

-0.07759 

22.000 

-0.22864 

A  l\r\rim 

VtUOilCHJ 

-8,88464 

23.000 

-0.53840 

0.03140 

-0.07090 

24.000 

-0.70549 

0.02258 

-0.04018 

Figure  8. 

Sample  listing 

of  tide  data  file 
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as  the  capability  of  computing  equilibrium  arguments  for  tidal  hindcasts 
or  predictions,  is  available  to  users  with  mainframe  or  PC  capabilities. 

The  intended  purpose  of  the  database  is  to  generate  time-independent 
realistic  tidal  data  for  a  specific  location  for  use  as  boundary  conditions 
for  the  long-term  fate  model  LTFATE.  However,  tidal  elevation  and  cur¬ 
rent  hindcast  or  prediction  capabilities  are  available  as  an  alternative  to 
published  tide  tables  that  give  only  high  and  low  tide  predictions.  The 
DRP  database  therefore  represents  an  improvement  to  the  tide  tables  be¬ 
cause  it  provides  a  capability  of  generating  a  continuous  tidal  signal  for 
any  time  period  at  any  location  in  the  computational  domain  shown  in 
Figure  1. 


References 


Borgman,  L.  E.,  and  Scheffner,  N.  W.  1991.  "The  Simulation  of  Time  Se¬ 
quences  of  Wave  Height,  Period,  and  Direction,"  Technical  Report  DRP-91- 
2,  U.S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 

Dronkers,  J.  J.  1964.  Tidal  Computatiotis  in  Rivers  and  Coastal  Waters,  North- 
Holland  Publishing  Company,  Amsterdam,  Holland. 


Luettich,  R.  A.,  Westerink,  J.  J.,  and  Scheffner,  N.  W.  1992.  "ADCIRC:  An 
Advanced  Three-Dimensional  Circulation  Model  for  Shelves,  Coasts,  and 
Estuaries;  Report  1,  Theory  and  Methodology  of  ADCIRC-2DDI  and  AD- 
CIRC-3DL,"  Technical  Report  DRP-92-6,  U.S.  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  MS. 


1QOO 
I 

Tn  f-iOrim  OffcKnrr 


"\AT 
r  it  . 


"a  a  ~ 

jc/iapcioiuiL  v 

n;  T  Ci  br- 


a.. 


^  1 :  ^ 


ictiyaib  ui  me  Kuiiiboldt  u&y, 

nl  O  1 tn  tari  m  t  Ht-e  horn  I  li  c'-nAnn  1  Cifo  ,r  Pn  1TYDT}  OO  1  T  T  C 

ah  Aruciun  wiioiuji  l  ivijyuoai  iDitc:,  iviidLtiiaiiuuuo  i  l^ivj  i ,  U.k^. 

Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 


Schureman,  P,  1940  (1958,  reprinted  with  corrections).  "Manual  of  Harmonic 
Analysis  and  Prediction  of  Tides,"  Special  Publication  98,  Coast  and  Geo¬ 
detic  Survey,  U.S.  Department  of  Commerce  (U.S.  Government  Printing  Of- 
fice,  Washington,  DC). 


Schwiderski,  E.  W.  1980.  "On  Charting  Global  Ocean  Tides,"  Reviews  in  Geo¬ 
physics  and  Space  Physics,  Vol  18,  pp  243-268. 

Westerink,  J.J.,  Luettich,  R.  A.,  and  Scheffner,  N.W.  1993.  "ADCIRC:  An 
Advanced  Three-Dimensional  Circulation  Model  for  Shelves,  Coasts,  and 
Estuaries;  Report  3,  Development  of  a  Tidal  Constituent  Database  for  the 
Western  North  Atlantic  and  Guif  of  Mexico,"  Technical  Report  DRP-92-6, 
U.S.  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 


'echnical  Note  OKP-1-13  (April  1994) 


13 


